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/^ ' Abstract 

Q> ■ The well-known Monte Carlo simulation packet JETSET is not built in order to 

G^ , describe few-body states (in particular at the few GeV level in e^e~ annihilation as in 
BEPC). In this note we will develop the formalism to use the basic Lund Model area 

CJ ' law directly for Monte Carlo simulations. 
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1 Introduction 

The Lund Fragmentation Model contains a set of simple assumptions. Based upon 
00 ■ them one obtains as a final result an area law for the production of a set of mesons 

Q ■ from a string-like force field. It is possible to reformulate this into an iterative 

cascade process in which one particle is produced at a time. This is essentially the 
way the model is implemented in the JETSET Monte Carlo packetW. 

In all computer programs it is necessary to make certain approximations, but 
Ph! the approximations in JETSET are hardly noticeable against the background noise 

Q^l signals as soon as sufficiently many particles are produced. For the few-body states 

(in particular at low energies) this is, however, no longer the case and in order to 
treat these situations it is necessary to use other means to implement the Lund 
Model as a Monte Carlo simulation program. In this note we will show how to 
make use of the basic area-law directly. We will be satisfied to treat two-body 
up to six-body states, which constitutes the overwhelming amount of the data 
obtained at the BEPC/BES accelerator in Beijing. 

We will start by briefiy reviewing some features of the string fragmentation 
scheme in the Lund Model. After that we will in the next section show how to 
implement the basic area law and present the Monte Carlo program LUARLW. 
Finally we will show some results, in particular point to some places where there 
are major deviations between the results from JETSET and LUARLW. 

Although the formulas in the Lund Model are derivedo, using (semi-) classical 
probabilities the final results can by comparison to different quantum mechanical 
processes be extended outside this framework'^. We will in this paper neglect 
all gluonic emissions and concentrate on the situation when the (color) force field 
from an original quark(go)-antiquark(go) pair, produced in e.g. an e"'"e~ annihila- 
tion event, decays into a set of final state hadrons. These are usually termed two-jet 
events. The transverse momentum of the final state particles will be treated in ac- 
cordance with a tunneling process, which leads to gaussian transverse fiuctuations, 
which are governed by the strength of the string constant k. 




Figure 1: String fragmentation and mesons production int — x space 

The color force fields are in the Lund Model modeled by the massless relativistic 
string with color 3(g) and 3(g) at the endpoints (the gluons {g, colorS) are treated 
as internal excitations on the string field). This means that there is a constant force 
field (k ~ 1 GeV/fm, corresponding to a linearly rising potential) spanned between 
the original pair. This pair is produced at the origin O, cf. Fig.|l|, and afterwards 
the (go^o) are moving apart along the x-axis (the longitudinal direction). The 
energy in the field can be used to produce new gg-pairs (new endpoints, thereby 
breaking the string). The production rate of such a pair (mass(es) fi, transverse 



momentum ±k_L, transverse mass(es) fj,± = J ji^ + k^ and with combined internal 
quantum numbers corresponding to the vacuum) is from quantum mechanical 
tunneling in a constant force field equal to 

The final state mesons in the Lund Model correspond to isolated string pieces 
containing a g from one breakup point (vertex) and a g from the adjacent vertex 
together with the produced transverse momentum and the field energy in between. 

In order to simplify the formulas and the pictures we will treat all gg-particles 
as massless and therefore as moving along light-cones (massive gg-particles move in 
this semi-classical scenario along Iwperbolas with these light-cones as asymptotes 
but the final results are the same'El). For the longitudinal dynamics in the Lund 
Model this is in accordance with quantum mechanics (in practice we only use 
dynamical features where position and rapidity of the particles are the relevant 
variables and this is allowed by the indeterminacy relationsQ). 

One necessary requirement is that to obtain real positive (transverse) masses 
all the vertices must have spacelike difference vectors. Together with Lorentz 
invariance this means that all the vertices in the production process must be 
treated in the same wayo. Another consequence is that it is always the slowest 
mesons that are firstly produced in any Lorentz frame (corresponding to the fact 
that time-ordering is frame dependent; this is also in accordance with the well- 
known Landau- Pomeranchuk formation time concept). Further each vertex has 



the property, cf. Fig. |1|, that it will divide the event into two jets; the mesons 
produced along the string field to the right and those produced to the left of the 
vertex. This feature implies that it is consistent to introduce a rank-ordering 
along the light-cones so that the first rank particle contains the original go (or go) 
together with the q {q) produced at the first vertex etc. Whether the process is 
ordered along one or another of the light-cones should be irrelevant. 

Based upon these assumptions it is possible to obtain a unique stochastical 
process for the decay ci'H. In particular this process will turn out to have simple 
factorization properties so that for every finite energy, rank-connected group of 
final state particles (which may possibly be part of the result from a large or even 
infinitely large energy reaction) there is a factor corresponding to the probability 
to produce the cluster with a certain energy and another corresponding to the 
probability that it decays into just this particular final state. We will from now 
on concentrate on this latter probability distribution. 

We find that the (non-normalized) probability for the cluster to decay into the 
particular channel with the n particles {pj} is given by 

n n 

dPint = lU Njd%]5{p'^ - m])5{Y.p, - Ptot) eM-bA). (2) 

3=1 3=1 

We note the appearance of the phase space for the final state particles multiplied 
by the exponential area [A is the shadow region in Fig. |l|) decay law. The 
quantity Ptot is the total energy momentum of the cluster so that Pf^^ = s, b 
is a fundamental color-dynamical parameter and Nj are normalization constants. 
Actually it is possible to generalize this result and include some special vertex 
factors but just as it is in connection with JETSET we have found no necessity to 
use such generalizations for the results in this paper. 

We will start the investigation assuming that there is a single particle species 
with mass m. Due to Lorentz invariance the integral of the distribution dPint can 
only depend upon s and this means that we can define the function i?„ for the 
multiplicity n by 

dPint{Pl,---,Pn) = Rn{s). (3) 

It is possible to subdivide these integrals into one part, that will depend solely 
upon the first-rank particle pi and a remainder part. The phase space element 
can be written as dpi5{p\ — m\) = dzi/zi with p+i = ziP+tot and p_i = m^/p+i. 
Similarly we find for the area 

A{lr--n) =ml/zi + A{2,---n). (4) 

Finally the energy momentum Pi, obtained after taking away the first particle, 
fulfills 

(Pi)2 = si = {Ptot - Pif = (1 - zi){s - ml/zi). (5) 



Therefore we obtain an integral equation linking Rn to i?„_ 

dzi 



Rn 

■ ' ■ Zi 



/iV,^exp(-6m?M)i?„_i(si). (6) 

J Z\ 



If we define the function -R(s) = I]-R„(s) the Eq (^ will be valid also for R itself 

n 

(to be precise there is a lower cutoff in z to make it possible to produce more than 
one particle). The equation has an (asymptotic) solution i? oc s" if 

1= f Ni — {l-z)^ exp{-bml) = f f{z)dz. (7) 

We have neglected a finite energy correction term (1 — m\/ zs)"" in the asymptotic 
limit. This defines the fragmentation function / in the Lund Model and this is the 
way JETSET implements the model. 

To be more precise, JETSET contains the production of all known particles. 
It is assumed that the parameters h and a are fixed by experiments (ordinary 
default values in the LEP region are h ^ 0.6 GeV~^ and a ~ 0.3 — 0.5) and that 
Nj normalizes the fragmentation function for each given mass irij. The relative 
probability of the different particle flavors is decided by an external procedure. 
Thus the production of strange quark-pairs (ss) are suppressed compared to u 
and d (cf. Eq.(|lD, an ordinary default value for the suppression factor is 0.3). It 
is further assumed (and this can be verified on a qualitative level) that it should 
be more difficult for a quark to tunnel into a vector meson (tensor meson) state 
(in the string model essentially a one- dimensional bag) than into a pseudoscalar 
state. For the production of baryons there are several more parameters (and recent 
investigationsHl indicates that for a description of baryon production one needs to 
introduce special features). 

There is another possible interpretation which has been pursued by a group 
from UCLAB. In that case all the numbers 6, a and A^ are taken as experimentally 
determinable parameters and the size of the fragmentation function integral in Eq. 
(|^) is used to determine the relative weight for the production of different particles. 
They use some straightforward Clebsch-Gordan coefficients, take into account that 
strangeness and baryon number must be conserved (thus e.g. the production of 
pairs of strange mesons must be used to determine the relative weight) etc. and 
with very few parameters they seem to obtain a phenomenologically good fit to 
data. 

There is, however, in both these implementations the problem of how to end 
the cascade process. In JETSET (and UCLA is using a similar procedure) the 
process is pursued stochastically from both endpoints and continues until there 
is a remainder mass of a certain size. This remainder is then fragmented as a 
two-particle state. This "joining"-procedure is implemented with great care in 
JETSET but nevertheless provides some unwanted features in particular for few- 
body states and low energies which is the main subject of this note. In the next 
section we will present a different procedure. 



2 Formalism 

Our basic assumption will be that the area-law, supplemented with certain simple 
probabilities similar to JETSET, will be fundamental also for low-energy and few- 
body states. We will further use a gaussian transverse momentum spectrum and 
the procedure is outlined in Appendix 1. We will for definiteness order the particles 
in rank along the positive light-cone. 

A particularly simple case is the two-body decay because in this case energy- 
momentum conservation determines all the observables. We obtain immediately 
for the quantity R2 (using rrij, j = 1,2, for the transverse masses; note that 
transverse momentum conservation in this case implies that p^i = — p±2) 



R2 = N,N2 ^y^ / '^^ cosh{^X{s,mlml) ), (8) 

'A(s, ml, ml) 



where the A is defined as 

X(x, y, z) = x^ + y^ + z^ — 2xy — 2xz — 2yz. 

The two different terms correspond to the two different values that can be chosen 
for the positive light-cone fraction Zj = p^j/P+ 



s + ml- 


-ml±^X{s,ml,ml) 


s +ml - 


2s 
-mf =F ^JX{s,mj,ml) 



We note that the two solutions will fulfill the relation Zj~z~s = m?,. This evidently 
means that we obtain the corresponding negative light-cone fractions by exchang- 
ing the sign in front of the square root. Although the expression for R2 looks 
singular at the threshold \/s = mi +m,2i it is in practice finite and even vanishing 
if the transverse momentum fluctuations are taken into account. 

From the expression for R2 we may derive the expression for i?3 according to 
Eq. @ but it is actually more convenient to consider R^ as a density in the area- 
size ^3. Again using the positive light-cone fractions Zj {j = 1,2,3) we may write 
the energy-momentum conservation equations and the size of A as 



2.7^ = ^ 



2 2 

mf m, 



— + — (l-zi)+m^ = A (10) 

Zl Z2 



Performing the integral we obtain 



m, = UN r'--''^^''^ , („) 



i=i ^A 



where 



A = [is-A)iA-ml^-ml2-ml^) 

22 22 22i2^222 

- mj_^mj_2 - m^^m^j^ - m^^m^^\ - Asm^^m^^m^^. 

A closer investigation tells us that the area A in this three-body case will be limited 
to the region 
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and that the quantity i 


\ can be written as 








A = {At 


-A){A-A,)[{s- 


-A){A 


-u)- 


-V + ■ 


where 


we have introduced the notations 
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W), (13) 



In this way the limits on the area corresponds to two of the zeros of the quantity A. 
It is rather easy to see that the particle configurations corresponding to ^ = A^ 
corresponds to the cases when the vertices (the breakup points) both lie on the 
same hyperbola. This means in particular that all the possible areas for the three- 
body decay are in between the two hyperbolas. To see that it is useful to note 
that the square root in the definition of Af can be factorized as 



(s — UY — A{y + vW) = [\/s — mi — 1112 — ma] [\fs + mi + 1712 — ms] 

X [v^s + mi -m2 + m32][v/s - mi + m2 + m3](14) 

and that it is a direct generalization of the corresponding expression for the two- 
body distribution 

A(s,m^,m2) = [v/s-mi-m2][v^ + mi + m2][v^ + mi-m2][v^-mi+m2]. (15) 

There are two sets of solutions for the light-cone fractions for given masses, energy 
and area size A (just as for the two body case). They can be written as 

. = -B±V^ .... 
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Figure 2: The vertex V divides the system into two subsystems of "right-m,overs" (mesons 
1, ■ ■ ■ , ni) and "left-movers" (mesons ni + 1, ■ ■ ■ , n). 



where 



B = A(A — m\^ — 771^2 ~ ''^±3) ~ "^(^ + ^"li ~ ""^±2 ~ ""^±3) 

+ m\^m\2 + rn^^m^r^ + m\r^m\^, (17) 



the corresponding values of Z2 and z^ can be read out from Eqs. (|10[). 

For the general n-body case we have to generate the kinenaatical variables by 
means of an (n — 2)-dimensional density but as we are going to limit ourselves to 
at most six particles we will make use of a different method. We will then for the 
four-body, five-body and six-body cases subdivide the system into two parts, each 
containing ui and n2 particles (^1,77,2 = 2 or 3 and n = ui + 722) . We may then, 
using the simple factorization properties of the model, apply the analytical results 
obtained above for the two parts. 

In more detail the subsystem closest to the positive light-cone (the g-end of the 
string) will be characterized by the positive light-cone fraction Z\_ = Zi and the 

(9) 

other subsystem by the (negative) light-cone fraction Z_ = Z2, cf. Fig.|^. Then 
the squred invariant masses of the two systems are 

si = Zi{l - Z-2)s and S2 = ^2(1 - ^i)s, (18) 

the area can similarly be subdivided into 

A = Ai + A2 + r with r = (1-Zi)(l-Z2)s. (19) 

In this way we obtain immediately for dRn the result 

dRn = ^^[exp{—bT^) + exp{—hV^)]dRnidRn2 (20) 

^JX{s, 81,82) 



in terms of the two possible F values 



r± = [s - si + S2 ± ^JX{s,Sl,S2) ] [s - S2 + si =F \JX{s,Si,S2) ]/4s. (21) 

The corresponding values for the Zj are (cf. the remarks after Eq. (P)) 



Zi = [s - S2 + Si±JX{s,Si,S2)]/2s, 



Z^ = [s-si + S2T^/Xis,Si,S2)]/2s. (22) 

We note that the momentum and energies of the two subsystems as well as their 
c.m.s-rapidities are easily expressed in these variables. Consequently the whole 
events can be determined in a straightforward way as soon as the (transverse) 
masses and the c.m.s energy is given. It is, however, necessary to introduce a 
method to chose the flavors of the different particles. 

We will use a method "in-between" the JETSET and the UCLA implementa- 
tions and in particular normalize our results at the c.m.s-energy ^/s = 4 GeV to 
the JETSET ones. Thus we will assume that the probability to obtain a partic- 
ular state with the n mesons labeled {rrij} (using the masses to characterize the 
particles) is given by 

P„(mi, ■ ■ ■ , m„; s) = C{VPS){SUD)N^R^{mu ■ ■ ■ , m„; s). (23) 

It turns out that the functions i?„ are the same independent of the rank-ordering 
of the particles. The quantity C is then introduced as a combinatorial number 
(stemming from the quark contents of the possible strings-there are in general more 
than one such ordered "quark-string" possible to make the particular mesons and 
even sometimes the possibility to make the same meson(s) at different places in the 
same string). The quantities (VPS) are the vector to pseudoscalar rate (we neglect 
the tensor mesons, although they can be introduced by straightforward means). 
The parameter [SUD) is the strange to up and down {qq) probability. It comes 
without saying that both (VPS) and [SUD) depends upon the number of choices 
to be made. The possibility to produce charm in the fragmentation process is-just 
as in JETSET- neglected in accordance with the results from Eq. (|T]). Finally A^ 
is a parameter to obtain the right multiplicity distribution. To obtain the right 
multiplicity distribution we will sum over all the different orderings and strings in 
Eq. (|23D . In this way we find a default value for the parameter A^ by comparison 
to JETSET. 

The procedure to produce a state is then that we firstly chose the multiplicity 
from the sum over all the Rn and then the particular state by a stochastical choice 
among the possible n-particle channels. After that we chose the particular ordering 
of the particles stochastically from the ways they can be produced from the possible 
quark-strings. Finally to obtain the particle distributions we just differentiate i?„ 
in the ways we have described it above. 

In this way we have defined a simple few-body low-energy version of the Lund 
Model (although an incorporation of baryons into the scheme may have some 



influence upon the results). It is implemented in the Monte Carlo simulation 
program LUARLWl^. A list of the default parameters of the model is given in 
Appendix 2. 

3 Some Results of the Model 

We will in this section exhibit a few results from the model implemented in the 
simulation program LUARLW, in particular discuss some configurations where the 
use of JETSET may be misleading in the few-body, low-energy cases. 

We firstly note that in a model of iterative character it is in general difficult to 
obtain the right behavior for small relative rapidities between two (rank-connected) 
particles. To see that we consider such a two-particle state, that is a part of 
a larger system, and assume that the pair has a total (transverse) mass-square 
Si2 = mf + ml + 2mim2 cosh.{6y) with rrij the individual (transverse) masses and 
6y the rapidity difference. This means that the distribution in 6y is 

dP dP ^ . wr N 

-TT- = - — 2mim2 smhfdj/J (24) 

ddy as 12 

Therefore unless the distribution in Si2 behaves (for values of ^/Sl2 close to mi+m2, 



i.e. for small 6y) as l/y A(si2,mf,m2) = l/2mim2smh.{6y) then dP/d6y will 
vanish. From the results discussed above this is in general taken into account 
by the generation mechanism in LUARLW but it is not so if we consider the 
corresponding situation in JETSET (cf. Fig. ^ and Fig.^. On the other hand it 
is difficult to observe this effect in a multi-particle environment so it is only for the 
few-body case it is of interest, especially if the phase space is small (low energies). 

There is further the "joining mechanism" in JETSET and related simulation 
programs. One place where this is noticeable is if we consider the rapidity ordering 
of a particular set of particles as they are generated in LUARLW and in JETSET. 
As examples we have chosen the four-body states p+p^vr+Tr^ and K^K^n^ii^ . 
In Table 1 we present the probabilities for the different possible orderings in two 
systems. 

We have only pinpointed a few particular places where it may be possible to 
see the precise workings of the area-law as compared to an iterative production 
mechanism. We find several places where there are clear differences but our major 
finding is that it quite surprising how well the basic JETSET scenarium really 
works at such low energies and multiplicities. 



Table 1: The probabilities (%) for the different possible rapidity order in JETSET and LU- 
ARLW for p+p^TT+TT^ and K^K^tx^tx^ final states. 



final state 


P+P TT^TT 


K+K-n+n- 


rapidity order 


JETSET LUARLW 


JETSET LUARLW 


yi> y2> y3> Vi 


13.15 


19.78 


18.55 


23.10 


yi>y2>yA> ys 


12.85 


12.39 


14.76 


12.85 


1/2 > 1/1 > 1/3 > 1/4 


12.60 


11.94 


14.70 


12.42 


1/2 > 1/1 > 1/4 > 1/3 


12.39 8.79 


8.47 


8.86 


1/1 > 1/3 > 1/2 > 1/4 


11.46 


12.43 


12.90 


12.82 


1/1 > 1/3 > 1/4 > 1/2 


5.98 


5.89 


3.72 


4.98 


1/3 > 1/1 > 1/2 > 1/4 


5.96 


5.96 


3.77 


4.83 


1/2 > 1/3 > 1/1 > 1/4 


5.29 


4.72 


6.67 


4.46 


1/1 > 1/4 > 1/2 > 1/3 


5.27 


5.11 


6.61 


5.33 


1/3 > 1/1 > 1/4 > 1/2 


3.48 


1.49 


0.41 


0.09 













Appendix 1 

The transverse momentum of the final state hadrons 

In the Lund model, a quantum mechanical tunneling effect has been used for the 
generation of quark-antiquark pairs g^gj. In JETSET, the p_L of a meson Qi-iQi is 
given by the vector sum of the transverse momenta of the gj_i and g^ constituents. 

In our scheme we produce the distributions of the particles directly from the 
area-law and then the transverse momenta {p±j} of the n particles must be de- 
termined at the same time as the longitudinal momenta {plj}- It is necessary to 
conserve the total transverse momentum and it is then necessary to include that 
in the generation. We have used two different methods which are described below. 

A Scheme I 



In this scheme we generate the transverse momentum of the particles directly, 
although one at the time. It is necessary to keep to total transverse momentum 
conservation and then it is necessary to modify the gaussian distribution. To that 
end we define the distribution 



_p{")| 



P±l? ■ ■ ■ ,P±n) 



'^(Epj 






rriy, 



pi, 



JJ exp{- 



i=i 



PJj. 
'2^2 



(25) 



where p±j is a two-dimensional vector, and a is the variance of the distribution. 

To obtain the distribution in p_Li, we must integrate the distribution F^"' over 
all the other vectors. We need to calculate the following integrals when we use the 



condition- density scheme to sample {p±j}. The distribution of p^i regardless of 
values of other transverse momentums reads, 

fl"'{Pll) ^ I '1^2- ■ ■ dWnF'-"Hpil,- ■ ■ ,Pln) (26) 

Under the condition of p^i is fixed, the distribution for px2 is 

/i"'(p±2|p±i) = /d'p±3 ■ ■ ■ d=Px,.-F"''(P±i, ■ ■ ■ , P±») (28) 

= N.-eM-i^^^^ii^^^^l^^^^]. (29) 

Similarly if all the vectors up to the j — 1th have been determined then the distri- 
bution for the jth {j < n — 1) particle is 

fj'^\p±j\P±l, ■■■, P±j-l) = J d^P±j+l ■ ■ ■ dp_LnF^''\p±l, ■ ■ ■ , p±n) 

(n-1 + 1) (P^7+'^'p±,/(n-J + l))2 

The final vector p_L„ is evidently determined by energy-momentum conservation. 
The effective variance of transverse momentum for jth particle is 



-f '" = V^r^TTT"' (^ = i.-- ■."-!)• (31) 

We see that ctj^Ii > (^j \ so the variances will be different for the particle with 
different order- number. This is because we have aleady determined a set of {p±}. 
If we instead ask for the inclusive distribution in px (i.e. the gaussian width 
independent of the rest of the particles), then it is the same for all values of j 



B Scheme II 

The method in scheme I is not the most common approach. In this case the 
transverse momentum conservation is fulfilled by construction and the particles 
obtain their transverse momenta from the constituents. At each production point 
the (gg)-pair is given ±q and the particle momenta are then 

p±i = qi, ■ ■ ■ , p±j = qj - qj-i, ■ ■ ■ , p±n = -q™- (32) 

There is a possibility (and at least for small mass-particles like pions this seems 
to be experimentally the case) that there is a correlation in the generation of the 



transverse momentums between adjacent production points, The most general way 
to introduce such a correlation in a forward-backward symmetric shape is, M 



n Y 



2a^ 



2 , (q2-P2qi) , , (qn-Pnqn-l) 

qi + — \ — 12 — + ■ ■ ■ + — \ — 2 — 



1-P^ 1-p^ 

(33) 
To see the symmetry we note that it can just as well be written in the following 
shape 

" 1 

/(")(qi, ■ . . , q„) = [] Si^jex-p - —^ q? + q^ + I] ^j(q? + q?_i - Se^qj ■ q,-_i) 



with 



4^2 



A, 



(34) 



(1 + p]) 


il-pj) 

2p, 



In the original paper ^, the correlation was phenomenologically taken as a func- 
tion of the mass of the particle. In general pj and therefore also ej are small 
numbers and can be neglected as they are in JETSET default. 



Appendix 2 

The parameters used in LUARLW 

In our notations, |Mp = exp{—bAn) is understood as the squared matrix ele- 
ment. The exclusive probabilites of n-body final state are given by 

Pn = J2[d^\M\\ (36) 

the following factors are considered in calculation. 

• VPS: A particle with spin J has 2 J + 1 spin projections, which means that the 
ratio of the vector mesons to pseudoscalar mesons should be 3 : 1. But e"'"e~ anni- 
hilation experiments show this ratio is smaller than the predicted value. There is a 
dynamical reason for the vector meson suppression e|l . Here, we treat the product 
of {2J + 1) and vector meson suppression factor as free parameters, and they may 
take different values for p/ir, K*/K and D*/D ratio. 

• C: In the quark model of hadron, a specified final state may produce in the frag- 
mentation of several different strings. We have counted the corresponding numbers 
C for all channels. 

• SUD: Lund model invokes the idea of quantum mechanical tunnelling to gen- 
erate the qiQi, which gives the relative probability for uu : dd : ss ^ 1 : 1 : Vg. We 
take r, ~ 0.3 same as in JETSET. 



• A^: As a default value, A^ is determined by requesting LUARLW to give out the 
same multiplicity distributions as in JETSET, 



N 



\ 



Y] pJETSET 



nfixed 



E CiVPS)iSUD)Rr^' 

nfixed 



(37) 



Table 2 The values of parameters used in Eq. (^3]) which are fixed 
at y/s = 4 GeV by comparison to JETSET. 



b 


VPS 


SUD 


N 




p/vr K*/K D*/D 




n = 3 


n = 4 


n = 5 


n = 6 


0.58 


1.07 1.15 


0.33 


0.49 


0.38 


0.26 


0.19 
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Figure 3: The light-cone component and rapidity distributions of a pirirp final state at 
s/s = 4 GeV (the ordering corresponds to rank- ordering). The distributions for the corre- 
sponding rank particle in JETSET (the net-like regions) and LUARLW (the solid lines) show 
significant differences, while the total inclusive distributions agree. The rapidity- difference in 
JETSET shows the gap near Ayij ~ (the so-called "disease" distribution). Figure (a) is 
non-normalized. 




Figure 4: The distributions of rapidity and rapidity differences for some special final states at 
/i = 4 GeV predicted by JETSET (the net-like regions) and LUARLW (the solid lines). 



